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Abstract 

The breakdown of conductance quantization in a quantum point contact 
in the presence of random long-range impurity potential is discussed. It is 
shown that in the linear response regime a decisive role is played by the 
indirect backscattering mechanism via quasilocalized states at the Fermi 
level; this can provide much higher backscattering rate than any direct 
backscattering process. For the realistic contact lengths (< 2000nm) the 
scattering processes prove to be independent, in spite of coherence of the 
electron wave. The distribution function of conductance fluctuations is 
obtained by direct numerical calculations as well as estimated within an 
analytical model for the first time. It is shown to be a generalized Poisson 
distribution. Estimates are made for quantum point contact performance at 
different choice of parameters. In particular, it is the better the larger the 
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intermode distance is compared to the amplitude of the random impurity 
potential. 

PACS: 72.10.Bg, 72.20.Dp, 73.20.Dx 
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I. INTRODUCTION 



The effect of conductance 2e 2 //i-quantization in quantum point contacts (QPC, elec- 
trostatically defined junction in a two-dimensional electron gas in a high mobility GaAs- 
AlGaAs heterostructure), first observed in 1988 [|l|, is still challenging both theoreticians 
ant experimentalists B. An intriguing feature of the effect is that it is not destroyed by 
elastic scattering, be it impurity scattering or scattering by a contact boundary 

The insensitivity of quantization to the latter process was explained within the frame- 
work of the adiabatic approximation [HI. Breakdown of quantization was shown to be 
insignificant as long as the contact shape (i.e., the confining potential induced by the gate 
electrode) is smooth on the scale of Fermi wavelength, A^. The last condition is likely to 
be true, since Xp ~ 40nm in a clean contact is at least an order of magnitude less than 
the characteristic scale of the gate potential variation [|J. 

The adiabatic approximation yields that different transverse modes pass the QPC 
independently ("no-mode- mixing" regime). This provides a sufficient, but not necessary, 
condition for the conductance quantization in a QPC: as was pointed in @||, the necessary 
condition is the absence of backscattering. The direct numerical calculations || show 
that conductance is quantized even if the intermode mixing is significant. The sum rule 
suggested there to explain this result follows from the unitarity of the scattering matrix 
of the system [f?j]. 

The impurity potential is more likely to produce backscattering, and thus to break the 
conductance quantization, since it is less smooth than the gate-induced potential. Many 
theoretical papers dealing with impurity scattering in QPC ||,[|] used the Anderson model 
with on-site disorder. This approach gives a qualitative understanding of the process, but 
its results cannot be directly applied to real GaAs planar structures, where the Coulomb 
impurity potential certainly is not sharp on the scale of A^ [fi|,|10|. In fact, the numerical 
calculations made for this realistic case [|J show that its scale of variation is intermediate 
between A^ and the size of the QPC itself. 

The effect of such a slowly varying random potential ("soft disorder") on the per- 
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formance of a QPC was investigated in f4,10[|. Its main feature is the negligible rate of 



scattering processes with large momentum transfer. Glazman and Jonson |TD] built a 
theory of the electrical conductivity through the QPC with soft disorder and showed that 
the direct backs cattering occurs mainly in the largest-number transverse mode (i.e., one 
with the minimal longitudinal momentum). On another hand, Laughton et al. || demon- 
strated that the long range impurity potentials may produce quasilocalized states inside 



the channel [|ll[ . Once appeared, such a state provides an effective indirect backscattering 
mechanism. Indeed, the quasilocalized state contains both "forward" and "backward" 
waves, so that the transition between the electronic states propagating in opposite di- 
rections via the quasilocalized state does not demand a significant momentum transfer. 
Thus, in an appropriate impurity configuration, indirect backscattering (formally due to 
the second order process) can have much larger cross section than the (first order) direct 
backscattering. The probability to find such a configuration of impurities, of course, must 



be investigated as well. Recently Gurvitz and Levinson |T2[ built a theory of resonant 
reflection and transmission in a QPC with a single attracting impurity at the contact's 
bottleneck. They have investigated the case when the first mode is open and a second 
one is being opened, thus creating a single resonant level at the top of the potential well. 
They obtained a Breit-Wigner type expression for the correction to the conductance at 
the transition region from the first to second conductance step (as function of the Fermi 
energy); its sign depends on whether there takes place tunneling into the resonant state 
(conductance enhancement) or scattering into it (suppression). 

In this paper we study numerically and analytically the indirect backscattering effects 
on the electrical conductivity of a QPC, in the presence of the screened Coulomb potential 
from randomly distributed charged impurities. 

In Sec.|IT] we find the corresponding correction to the current, supposing that there 
are quasilocalized states in the QPC, and that the impurity potential is soft, so that 
the tunneling processes to and from these states can be neglected, and their energetic 
spectrum is dense. We show that the correction, given by a Breit-Wigner type formula, 
suppresses the current. If the quasilocalized states exist close to the Fermi surface, this is 



really the leading contribution from the impurity scattering in a QPC. The result is valid 
for the bulk of the conductance step vs. gate voltage. 

In Sec.[TTl] we investigate the stochastic properties of arising conductance fluctuations, 
as a random variable determined by the specific impurity arrangement. The distribution 
function of the conductance fluctuations is shown to be of Poisson type and is explicitly 
determined by the tunable parameters of the QPC: its length and number of open modes 
(i.e., conductance). This allows us to make more detailed predictions about the QPC 
performance than mere knowledge of the fluctuation dispersion. 



In Sec.|V| we estimate the theory parameters based on the quasiclassical description 
of the random impurity potential. 

In Sec.[V| we present the numerical calculations that provide the basis for the present 
work. The impurity potential was calculated in a self-consistent way using the Thomas- 
Fermi approximation. The quantum contact was modelled by imposing upon it a parabolic 
confining potential. The conductance was calculated in the linear responce limit by a 
standard transfer matrix method for different contact lengths and for different realizations 
of the impurity potential. ( Different realizations of the random impurity potential in the 
QPC was performed by changing the position of the confining potential with respect to 
once obtained random potential profile.) 

An empirical distribution function of conductance fluctuations was thus obtained, and 



we compare the results of the numerical calculations and theoretical predictions in Sec.[V]. 
They are in very good agreement. This allows us to conclude that the leading mechanism 
of conductance quantization breakdown in QPCs is the indirect backscattering, and it 
allows us to find estimates of the QPC performance. 

II. INDIRECT BACKSCATTERING IN A QUANTUM POINT CONTACT 



In accordance with what is now the standard approach ||[10|, we start from the adi- 
abatic model of a QPC, where the transverse modes are well defined; the deviations will 
be regarded as a perturbation (see Fig.Q). The Hamiltonian of the electron contains the 
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following terms: 

H = -&(&- + ^) + u ' [x ' y) + v{x ' y) - (1) 

Here U g (x,y) is the confinement potential, shaping the junction in a two-dimensional 
electron gas (2DEG). It proves to be convenient to include the smooth component of 
the impurity potential as well. The rest of impurity potential, V(x,y), is regarded as a 
perturbation, e.g., leading to intermode mixing. It is also a soft potential, which does not 
lead to large momentum transfer. 

The potential U g (x,y) is supposed to be a slow function of the coordinates. Then the 
wave function of an electron with energy E inside the QPC can be expanded as follows 
I: 

y;E) = J2 Cm,ai) m ,a( x > V\ E )- ( 2 ) 
The adiabatic eigenfunctions are defined by 

ip m , a (x, y; E) = cp m (x, y)Xm,a( X 'i E )i ( 3 ) 

where the transverse eigenfunction corresponds to the m-th eigenvalue of the transverse 
Hamiltonian: 

#L0m(£,y) = \ ~ 2m*d 2 + U o( x iV)\ < l > riS?,y) = E m± (x)(p m (x,y), (4) 

and the longitudinal eigenfunction satisfies the equation 
f h 2 d 2 ) 

I ~2m*dx^ + Em± ^ [ Xm > a ( x ' E ^ = E Xm,a{x\ E). (5) 

We can distinguish different groups of electronic states in the QPC, according to the 
behaviour of Xm,a(%i E) in the presence of the effective one-dimensional potential E m ±(x) 
(see Fig.^). They are denoted by the index a = —2, — 1, ..2. 

The states of special interest for us are the quasilocalized ones (a = 0). They ap- 
pear if the effective potential E m j_(x) is a nonmonotonous function. Note that for each 
quasilocalized state in m-th mode, generally, there exist propagating states with the same 
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energy in some mode of less number (Fig|J). This is the reason why we these states 
are gnaszlocalized (due to smoothness of the potentials in the QPC we can safely neglect 
another reason for this, namely, the tunneling decay). 

The coexistence of propagating and quasilocalized states (P-states and Q-states) at 
the same energy is characteristic for QPCs and gives rise to indirect backscattering in the 
QPC. 

The correction to the current due to backscattering can be written in a standard way 

IB 

A/ = -HM J dE {n F {E - m) - n F {E - fi 2 )) J(E); (6) 



N N 



^) = EEM£)|. (7) 

m=l n=l 

In the linear response limit this yields the correction to the conductance, 

AG = -^J(E F ). (8) 

The summation in (|7|) is taken only over N propagating (open) modes in the quantum 
point contact, and r nm (E) is the probability amplitude of the electron incident from the 
right in mth P-mode (a = 1) to be scattered back, to the nth P-mode with (a = —1): 

r nm (E)5(E - E') = (^(E'), ■ (9) 

Here ty™\(x, y; E) is the scattered wave, corresponding to the unit wave with energy E, 
incident from the right in the mth mode, ip m) i(x, y; E). This function can be found from 
the Lippmann-Schwinger equation jL4j] (r = (x,y)): 

#™i(r;£) = Vm,i(r;£) + / dr>G R (r,r'; E)V(r')^ mA (r>; E). (10) 

In this equation we have introduced the exact retarded Green's function (retarded Green's 
function) in the presence of the perturbation V(x,y), for which we can write : 

G*(r, r'; E) = G*(r, r'; E) + J dr"G*(r, r"; E)V(r")G«(r, r'; E) + ... (11) 



The unperturbed retarded Green's function is given by 

G*(r, r'; E) = £ ^ig^ij s £ G *«( r , r ' ; £). (12) 

m,a,e t -t- «U Q 

The summation is taken over all the quantum numbers of electronic unperturbed eigen- 
functions. 

Inserting (|10 JTlf) into @ we find two first nonvanishing terms: 



r nm (E)5(E - E') = (ip n ^(E'),G^(E)V^ mA (E)) + {^ x {E'),G^{E)VG^{E)V^ mjl {E) 

We neglect the higher-order terms which is equivalent to an assumption that the scattering 
processes by different quasilocalized levels are independent. This assumption is valid 
only in the case L <C k r ,hoc, where k r ,hoc are the transport and localization lengths, 
respectively. The following results show that this condition holds in the case of realistic 
quantum contacts. 

According to our initial assumptions, the perturbation potential is soft. This means 
that, generally, its matrix elements between the P-states propagating in opposite direc- 
tions (i.e., with a = +1 and —1) are negligible, while the matrix elements between the 
P-states and Q-states (a = ±1 and 0) are nonzero. Since the retarded Green's function 
(|12|) contains the products of electronic eigenfunctions with the same indices a, Eq. (|l3D 
reduces to: 

r nm (E)6(E - E') = ty ni - 1 (E'),G*-'-\E)VG*'°{E)V'il; m , 1 (E)) . (14) 

The superscripts in the retarded Green's functions show which part of the expansion (|12"D 
we keep. 

The formula ( |T4"D shows that the backscattering from the incident (m,l)-state to the 
(n,-l)-state occurs through the set of quasilocalized states (described by the part of re- 
tarded Green's function, denoted by Gq'°(E)). 

After some standard transformations, we obtain that 

r m JE) = 2ME) E ( " ; ^ E] V k E °':] ( f§ C| V Im; l! E) ■ (15) 



Here we have introduced the one- dimensional electronic density of states (DOS) at infinity, 
v(E), to account for the propagating modes, and denoted the matrix elements of the 
perturbation potential by 



1; E\ V \q; 0; e) = J dr^^r; E)V(r)^ (r; e). (16) 



The summation is now taken only over the quasilocalized states; as such, they have 
discrete spectrum, so that e is a discrete variable; q is the number of transverse mode 
where the localized state appears. 

The equation (I5Q contains the unperturbed retarded Green's function of the electron 



in localized states, which thus have an infinite lifetime. The perturbation, mixing different 
states, makes the Q-states metastable, and instead of iO-term there appears the spectral 
function iT q (e, e)/2 [0 (the accompanying shift of the energy levels can be accounted for 
by changing the summation variable e). The spectral function is given by 

r 9 (e,e')=27r£ £ \(q; 0; e| V \n; a; e')\ 2 u{e); (17) 

n a =±l 

(of course, a Q-state can decay only into P-states). 

The correction to the current (§) is expressed through J(E) = Y, m ,n \ r nm(E)\ 2 . Note 
that the matrix elements corresponding to the Q-states from different subbands (q ^ q') or 
different impurity potential wells in the same subband (w ^ w', where w is the label of the 
potential well) enter the expression for r nm (E) with their phases, which are uncorrelated. 
This means that the main contribution to J{E) will be given by the diagonal terms in 
the corresponding sum, that is, as it is easily seen, 

J(E) ~ 4vrV(£) x 

v v E m ,n (n; -l;E\V \g; 0; e) (g; 0; e| V \m; 1; E) ■ (n; -l;E\V \q; 0; eQ (g; 0; e'\ V |m; 1; E) 
bb (E-e + iT q (e, e)/2) (E - e' - iT q (e>, e')/2) 

Here e, e' are different quasilocalized states in the same potential well of the same ID 

subband. 

This expression can be further simplified, if the distance between the energy levels in 
a potential well are large comparatively to their width: Ae we u ^> T. Then in the sum 
over e, e' in (0) the terms with e/e' are negligibly small, so that finally we obtain: 



2 2( ^ ^ \(n; -1; £| ^ |g; 0; e)| 2 Kg 0; e| V \m; 1; £)| 2 
1J " ( j t? (E-^ + (r„( ei c)/2)« 

= y g^gVg)! (1Q) 

The result has a typical Breit-Wigner form for a set of independent resonant levels 



14| , which is the case under the assumptions made. 

Though the energy spectrum of the Q-states in a single potential well is assumed to be 
rarified, in a long contact, where many potential wells with random parameters appear, 
and due to the fact that each Q-state has a finite width T, the spectrum of these states for 
the whole system is dense enough to be described in countinuous approximation. Then 
we can introduce the density of localized states in the qth mode, N 9 (e), and rewrite (|T9| ) 



as 
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In the linear response limit this directly gives the correction to the conductance: 

AG ~ "¥i £ **@F)r q (E Fi E F ). (21) 

We discussed only virtual (including resonant) scattering to and from the Q-states. In 
reality, when the finite driving voltage U is applied, there exist the processes of real elastic 
and inelastic (electron-electron or electron-phonon) scattering between Q- and P-states 
with energies in the interval eU around the Fermi energy. The contribution of these states 
to the scattering is proportional to eU. We can safely neglect their existence in the if 
the linear response limit eU — > 0, when the contribution from the Q-states outside the 
et/-band (in the interval of width ~ T(Ep, Ep)) is dominant. These kinetic processes, of 
course, must be taken into account if we would like to discuss the nonlinear response of a 
QPC. 



III. STATISTICS OF CONDUCTANCE FLUCTUATIONS 



Now let us discuss what the formula ( plf ) yields. The conductance quantization break- 
down can be characterized by the relative difference between the actual conductance 
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vs. gate voltage dependence, G(V g ), and the ideal one, averaged over the n-th step 
(Fig.|,inset) (reduced conductance deviation): 

Inserting here fl2"T|), we get: 

9 "=2?4 W j W V, (n + 1 )-nw (23) 

where the n-th mode opens at the gate voltage equal to tg(n). In the last equation we can 
instead of integrating over the gate voltage integrate over the subband bottom energies, 
E±(x), which are determined by V g . The intermode distance, AE±(x), weakly depends 
on the gate voltage, the main effect of which was directly shown to be an upward shift 



of the potential in the constriction [nj. In the realistic model of a parabolic confining 
potential it also does not depend on the mode number, and in the adiabatic limit is almost 
coordinate independent. Its magnitude is of the order of E F /N max , where N max is the 
number of transverse modes in the wide part of the constriction (and equals approximately 
to the number of conductance 2e 2 //i-steps which can be observed in given QPC). This 
enables us to simplify the above equation: 

^^r^w^M. (24) 



The functions under the integral are implicit functions of the bottom position of the 
corresponding mode. 

If localized states have an almost continuous spectrum (which is consistent with our 
initial assumption of the impurity potential being smooth), the main dependence on the 
bottom energy of the mode enters the DOS, not the level width, which is a smooth function 
of energy in the scale of AE±. Therefore we can take it from under the integration. Then 
we obtain the following formula: 

9 «™2' AE L L q± ^ ? (25) 
Here we explicitly show that, as follows from the definition (|17D, T q (E F , E F ) depends on 
the number of propagating states in the contact, that is, on the step number n. 
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The integral over energy in fl25|) has a sense of number of localized states in all ID 

subbands in the contact, which pass the Fermi level as we sweep V g across the n-th 

conductance step, Nq. We will denote its average over possible realizations of the random 

potential by kL[L is the length of the contact): 

(No) 1 rE q ±+AE ± / \ 

I dE ± (Y,X q (E F )), (26) 



L JE ql 



or 

rE q± +AE A 
K = \ 



dE ± (^2M q (E F ) ) . (27) 



Here (M q (E)) is the ensemble average of DOS per unit length in the g-th ID subband over 
possible realizations of the random impurity potential, so that the ensemble average of g n 



is 



(9n) = ~ ■ ^^ kL. (28) 

Now the expression (|25|) acquires a clear physical meaning. The correction to the 
contact conductance due to indirect backscattering is a sum of contributions from inde- 
pendent Q-states (in each ID subband): 

g n = N Q Ag n . (29) 

Each contribution is approximately the same and is proportional to n (see (|1 



7i(T(E F ,E F ;n)) 
n = 2AE { = in - (30) 



The density of states is the self-averaging quantity ||15|| , so that in the limit of infinitely 
long quantum contact we would get : 

lim ^ = (M q (E)), (31) 

Nn 

lim = k. (32) 

L~^oo E 

In the contact of finite length the number of appropriate Q-states deviates from its 
ensemble average, (Nq) = kL, and thus the conductance of the QPC fluctuates from 
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realization to realization. In order to predict the performance of a single QPC, we need the 
distribution function of these fluctuations, dependent on such parameters of the contact 
as its length L and number of open modes n. 

On a length scale large compared to the size of the localized state (correlation radius 
of the impurity potential, of order 100 nm), we can regard these states as independently 
randomly distributed along the channel, with occurence per unit length k. The probability 
to find p such states along the total length L of the contact is then given by the Poisson 



formula 18 



U{N) = (kL) n e~ KL /N\ . (33) 

Since each state gives the same contribution to the conductance, Ag n = jn, then the 
probability density of the conductance deviation on the n-th step 

n9n = N ln) = LJ__ . (34) 

The reduced conductance deviation acquires the discrete set of values, g n = 
771, 27n, 3771..., depending on the number of relevant Q-states in the contact. 

The shortcoming of the above formula (discrete set of values of g n ) follows from our 
simplifying assumption that each Q-state has the same width (T(Ef, Ep)). In reality, this 
quantity also fluctuates, so that the variable g n is rather a continuous one. Nevertheless, 
as the numerical calculations show, this expression provides a fairly good description of 
the performance of a QPC (see Sec.[VI|). 



IV. ESTIMATES OF THE THEORY PARAMETERS 

We see that the statistics of the conductance fluctuations due to indirect backscattering 
is determined by the following parameters (see Eq. (|3^) ) : effective length of the QPC, L, 
number of open modes, n, average ID density of Q-states at the Fermi energy per unit 
length, k, and average contribution to the conductance from each Q-state, Ag. Since 
the latter quantity is proportional to the average scattering rate from these states to the 
propagating ones (see ([17D ), 

13 



it is more convenient to use the set of four independent parameters, L, n, k and 7. 

The first two are tunable parameters, which can be changed by changing V g and/or 
gate configuration, while the other two are essentially determined by the properties of 
given GaAs structure (density and charge of impurities, spacer thickness etc.), or, in 
short, by the impurity potential in the system. We need some numerical estimates for 
these intrinsic parameters, 7 and k. They can be obtained from such characteristics of 
the impurity potential as its correlation length, ly, and dispersion, a 2 = {V 2 mp ), which 
are both contained in the correlation function K(x) or spatial spectral density S(k) (see 
0): 

K(x) = (V imp (0)V imp (x)); 5(k) = J ^xK(x)e- ik ' x (36) 

(We put to zero the average value of the impurity potential.) 

First we estimate 7. Since (T{E F , E F )) ~ 2vr • n ■ 2v\(Q\V\P)\ 2 (see (0)), then 

7 ^ 2ix 2 v\(Q\V\P)\>/AE±. (37) 

Here |(Q|^|-P)| 2 is the average square modulus of the impurity potential between quasilo- 
calized and propagating states. Evidently, it is of the order of S(AkpQ)a 2 /l , where Ak P g 
is the difference of longitudinal wave vectors in P- and Q- states, a is the reduced matrix 
element of impurity potential between different transverse modes (a ~ 0.1 according to 
our numerical calculations, see Sec.0), and l ~ ly is the average longitudinal size of 
the quasilocalized state (the latter factor appears because the wave functions of Q-states 
are normalized to unit probability, while the P-states are normalized to unit probability 
current along the axis of the QPC). In accordance with our basic assumptions, for the 
relevant (Q,P)-matrix elements of the impurity potential the longitudinal wavevectors are 
almost the same: Akpq <C kp. Therefore we can take 

\{Q\V\P)\ 2 ~ S{0)e~ 2Akp * lv a 2 /l v = S{0)e~^ ^~a 2 /l v ~ a 2 l v a 2 . (38) 

This gives the following estimate for 7: 
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We use here the formula for ID DOS of propagating states, v = l/nhvF- 

If we insert into fl39|) the values consistent with our numerical calculations, a = 0.1, 
X F = 40nm, l v = lOOnm, a 2 = 0.1E F , AE ± = 0.3E F} then we obtain 7 < 0.15. This 
is a proper order of magnitude estimate, though higher than the value that follows from 
the numerical calculations (see Sec.[VT|). The reason is, that putting Ak P Q = in flSSp we 
overestimated 7. To eliminate this discrepancy, it is sufficient to take Akpq/kp — 7%, 
which is quite reasonable. 

One useful remark: In the superposition approximation (see, e.g., fll6| ) we can present 
the impurity potential as a sum of potentials of identical potential centres, v(x), with 
density n imp randomly distributed along the conducting plane. Then by the Carson 
theorem we get |L8| S(k) = ni mp \v(k)\ 2 . Therefore the parameter 7 is proportional to the 
impurity density. 

In order to estimate the other intrinsic parameter, k, we use the quasiclassical, or 
Thomas-Fermi, approximation for the average density of localized states in the ID sub- 
band, (J\f(E)) (see flEfl), leading to the following approximate expression: 



>/™ am t \ ( ( E i±( x ) ~ E )*\ ^ (E ql _(x)-E 



(Af q (E, x)) = -*-—6(E q± (x) - E) exp - v '- D_ v 



2 



a 



2- 



^ 6(E q± (x) - E) VTE _ f (E q± (x)-E) 
~ y/E q± (x)-E " I )■ m 

D -1/2(2) is the parabolic cylinder function. As a matter of fact, the validity limits of the 

latter expression are determined not by the quasiclassicity conditions, but by much looser 

ones 
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A F 



E-E q± (x)\^a,E F --f. (41) 



lv 



Substituting it into the definition of k (1271), we find 



^ +AE± , 0(E - Ep) sjm {E-E f f 



K ~ / dE - =^-^e 



^E - E F irh 



q 

00 dE . „ 1 Jm ( e - e f? 



AE ± =-^r e • ( 42 ) 



Ef+aei Ai?_L \JE — Ep irh 
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The latter formula can be identically rewritten as follows: 



K , 2 Wi,WY& i (43) 



4' 2 V a ) J \ E F X F 
Here T(a, z) = dy y 1 ~ a e~ y is the incomplete gamma function. 



With the same choice of parameter values as above, we find that k ~ V0.66 
r(l/4, l/2)/\ F = 0.269/A F « 0.006 nm" 1 . 

V. NUMERICAL CALCULATIONS 

The basis for our hypothesis about the statistics of conductance deviations from ideal 
quantization is the numerical calculations of QPC conductance for different realizations 
of soft random impurity potential, at different contact lengths and for different number 
of open modes. 

The model used for simulation of a heterostructure device is very similar to the one con- 
sidered in [|J. Donors are assumed to be fully ionized and distributed randomly through 
the donor layer. We restricted donors to a plane which should be considered as the middle 
of the donor layer. In our calculations all the donors at the n- type Al^Gai-^As- doped 
layer have the same height /i=30 nm above the 2DEG. The electrons are treated as a two- 
dimensional layer of 10 nm thickness, which is much smaller than other relevant length 
scales. The positions of the impurities were generated by the uniform random number 
generator at the square 1290x1290 nm 2 . We assumed periodic boundary conditions, i.e. 
this square was continued periodically in all directions in 2D plane. In order to avoid the 
occasional appearance of too dense clusters of impurities at the donor plane the distribu- 
tion of impurities has been 'relaxed', namely, we did not allow the appearance of more 
than 5 impurity atoms at the area 1/n, where n is the 2D concentration of impurities. 
For the realistic concentration of impurities, n = 10 12 cm" 2 one must discard not more 
than 5 -f- 10 % of impurities. The potential of a single impurity was taken in the following 
form: 

v{r - n) = . 6 =, (44) 
e v /(r-r l ) 2 + / i 2 
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where e (« 13) is the static dielectric function of GaAs , r is the position vector in the 2D 
plane, and rj is the coordinate of i-th impurity. Since we did not want to be too specific 
in modelling any definite heterostructure, we did not take into account the image term 
considered in ||. 

The potentials of the impurities were summed directly. The numerical calculations 
show that the summation of the potentials of impurities situated within the radius |r — 
ft I < r max = 10 h gives the resulting potential with a good precision. For the rest of the 
plane the summation could be substituted by integration over a plane with homogeneously 
distributed charge. This gives a constant term, which can be dropped since we are only 
interested in the fluctuations of the impurity potential around its average value. Therefore, 
we choose the average value of the summed potential as a point of reference for the energy, 
i.e. we take {V imp ) = 0. Further increase of r max does not give any significant change for 
the potential fluctuations defined in such a way. 

The amplitude of fluctuations of the unscreened impurity potential proves to be too 
large, for a realistic concentration of impurities n ~ 10 12 cm~ 2 being a few times greater 
than the Fermi energy (of the order of 10~ 2 eV). This means that the screening of this 
potential by the electrons in 2DEG should necessarily be taken into account. For a 
qualitative estimation of this effect we used the way proposed in [f|: the Thomas- Fermi 
approximation, which is applicable in the case of slowly varying (on the scale of A_p) 
impurity potential, i.e., our case. The density of electrons is then given by the local 
equation 



where fj, is the chemical potential, and Vtot is a sum of the unscreened impurity potential 
and the induced potential from the electrons in 2DEG, 



Vi n d is related to n in ^ by the Poisson equation. For a fixed chemical potential Eqs. fl45|) 
and fl4"6|) give a possibility to find the total potential V to t in a self-consistent manner. One 




(45) 




(46) 
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must keep in mind that the experimentally observed Fermi energy should be measured 
from the average total potential, namely, one more equation should be added to this 
self-consistent scheme, 

E F = /i-(V tot ). (47) 

Here (Vtot) is the average total potential, Ep is considered as a fixed parameter given 
by the experiment, \i is the "bare" chemical potential. It shows the degree of filling of 
the system by electrons and varies from iteration to iteration. Obviously, this system of 
equations has a single self-consistent solution. The numerical calculations show that for 
n ~ 10 12 cm -2 and Ep ~ 10~ 2 eV this solution is achieved when the system is almost 
completely filled by the electrons, i.e., for the most of the points in 2D plane ni n d is 
non-zero and the resulting potential fluctuations are smaller than (or of the order of) the 
Fermi energy, V to t — {Vtot) < Ep (see Figjj].) The characteristic scale of the change of the 
potential fluctuation is of the order of 100^-200 nm. 

The calculations of the potential were made with a crude mesh in 2D plane. To solve 
the scattering problem by the standard transfer matrix method, a finer mesh is desirable, 
and the potential was therefore interpolated by the cubic splines when necessary. 

In all these considerations we did not take into account the self-consistent effects of 
redistribution of electrons in the constriction due to the gate voltage. Nevertheless, it 
is clear that such effects could be very significant, since it could enhance the potential 
of a single impurity located in a bottle neck because the concentration of electrons in 
this area is lower than the average one, and screening is suppressed (impurity denudation 
effect). Such an effect increases the potential fluctuations at the bottle neck and can 
change significantly the whole behaviour of the conductance (locking the channel) 0. 
The question of how this effect of just a few impurities located at the narrowest place 
may show up in the statistics of conductance fluctuations is still open. We plan to discuss 
such a challenging subject elsewhere. 

In order to obtain the transmittance and reflectance matrices we solved the 2D 
Schrodinger equation rewriting it in matrix form by the standard transfer matrix method 
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with the proper boundary conditions. In order to model the gate voltage we used a 
parabolic potential, as in ||, 

U g {x, y) = hy^m^y/bix)} 2 , (48) 

where 77 = l/4fcp&(x), and 

b( x ) = - {boo - b ) sin 2 (7rx/L); -L/2 <x< L/2. (49) 

We chose the parameters bo=10 nm and &oo=170 nm, i.e., the maximum number of the 
modes passing through the constriction is eight, for the realistic Xf = 42 nm. 

As long as we are interested mostly in long constrictions (of length 600 nm and longer) 
this impurity scattering term will be dominating in comparison with effects of changes 
of the width b(x) along the constriction, a case thoroughly considered by Brataas and 
Chao || for short contacts. In the bulk of our calculations we neglect them. Because 
the most important region for the breakdown of quantization lies near the bottleneck, we 
multiply the impurity potential inside the contact by the factor sin 2 (7nr/L). This gives 
us the freedom not to care about changing the boundary conditions for any realization of 
the potential. The number of entering modes is always constant. 

The statistics of conductance deviations in quantum contacts can be obtained either 
(i) by generating the impurity potential (with the self-consistent procedure) for each re- 
alization of the constriction, or (ii) by moving the center of the contact (possibly also 
changing the direction of the axis of the contact) over the plane with the impurity poten- 
tial obtained in a self-consistent manner only once. The second way is surely less time 
consuming and, therefore, is the one we use. 

We performed numerical calculations for contacts of three lengths (600, 800 and 1200 
nm) for 625 different realizations of impurity potential, moving the center of the X-oriented 
contact along the square of 1290x1290 nm 2 , with the periodic boundary conditions de- 
scribed above. 

We also calculated the conductance of very long QPCs (up to 10000 nm) for several 
realizations. 
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VI. ANALYSIS OF THE RESULTS 



The results of the calculations are shown in FigsfJ,^] and |5|. Typical conductance vs. 
gate voltage curves are shown in Fig|| Curves A correspond to an impurity free contact. 
The center of the contact lies at the same place for all the curves B (at a minimum 
(valley) of the impurity potential) and C (at a maximum (hill)). Evidently the C-curves 
demonstrate better quantization than the B-curves, for the same contact length and step 
number. This agrees with our statement that the indirect backscattering is more effective 
mechanism of the quantization breakdown in QPCs. The quantization quality becomes 
poorer for larger contact length and higher step number. 

The conductance vs. gate voltage curves in a very long QPC (L=10000 nm) for three 
different realizations are shown in Fig.f|. We would like to note here that no trace of 
localization in the system is detected. On the other hand, the conductance plateaus 
(except the first one) are totally destroyed at this length. 

In order to obtain more convincing proof of our theory, the statistics of conductance 
fluctuations was investigated. The average conductance deviation from the ideal con- 
ductance step , g, was calculated in 625 realizations for each step (n = 2,3, ..,6) and 
each contact length (L = 600,800, 1200nm), and the empirical distribution function was 
obtained: 

Pe mp (g)Sg = N(p6g)/Q25, (50) 

where N(p6g) is the number of realizations for which p5g < g < (p + l)Sg. In the actual 
calculations 5g = 0.02. 

Since the contact has a finite length, the deviation g is finite even without impurities 



due to the geometric step smoothening [IS]; it is the larger the shorter the contact is. 



Therefore when calculating g, we have substracted from it the corresponding value of g 
for the impurity free case. 

The average deviation increases with the mode number and the contact length. The 
empirical distributions are essentially asymmetric. Some of them have a tail stretch- 
ing into the negative g region. This means that in some cases the impurity potential 
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can make the contact more transparent in comparison with the case without impurities. 
There may be two mechanisms for this transparency growth. First, there may occur 
tunneling through the localized states in the transition regime, which enhances the con- 
tact transparency fll2"| ; second, the random impurity potential may effectively lessen the 
aperture of the contact, thus making the conductance steps sharper and enlarging g |T9| . 



Both mechanisms are consistent with the fact that the negative tails disappear for larger 
contact length and/or step number, and both processes are likely to occur in different 
random potential arrangements. 

The values of g obtained from the numerical calculations are not quantized, as the 
expression ([34]) implies. This is, as stated above, due to the fact that the widths of 
different quasilocalized states are not exactly the same. Nevertheless, by interpolating 
this formula to noninteger values of g we achieve a strikingly good description of the 
numerical results: 

ULY^ n e~ KL , . 

™ = TO^TT)' (51) 

This result can be easily understood, if we take into account that the contribution of each 
single Q-state to the correction to the conductance is small enough ( a few percent of the 
total effect, see below), so that a comparatively large number of the indirect backscat- 
tering processes is necessary to obtain the average effect in the QPC. Then the discrete 
distribution (34[) almost coincides with its continuous interpolation (|5lD. This is the more 



true, the longer the contact is. Indeed, the fitting of the two curves is better for the larger 
values of L (see Fig.|). 

The solid lines in Fig.|5] are the least square fits of Eq.(|51|) to the empirical distribution. 
The dots on these lines indicate the points where g is an integer multiple of 777. (according 
to (|34|). We use the same value of the parameter k = (1/320) nm" 1 = 0.13/Ai? for all the 
curves, while the parameter 7 varied (the mean square error of the fitting was in the most 
cases practically the same as if we would vary them both independently). The choice of 
k agrees very well with the estimates of Sec.|lV| ( k ~ 0.269/Ai? ~ 0.006 nm -1 ) , up to a 
factor of 2. This shows that the effective contact length where the scattering occurs is of 
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the order of L/2, i.e., coincides with the length of its bottleneck part (see Fig.||). 

The values of 7 are shown in Table |. They change insignificantly with L and n, and 
are in a good agreement with the estimates of Sec.|V[ The theoretical curves provide very 
good description of both position and magnitude of the distribution peak, as well as of 
its large deviation tail. 

It is noteworthy that the width of the resonant peaks in G(V g ) curves (Fig.|]) agrees 
with our estimates of the resonance width T ~ '-/nAEx (AE± being of the order of the 
stepwidth). 

The knowledge of the distribution function of conductance fluctuations allows to obtain 
more accurate criteria of good performance of a QPC, e.g., to predict the probability for 
a QPC of a given size, realized in the GaAs structure with given properties, to have a 
certain number of well-defined conductance steps. 

The analysis of the theory parameters 7 and k shows that for different experimental 
situations they can be estimated as follows: 

7 ~ o.016-^; (52) 

K(nm -^_L.exp(l--^-(^-) j. (53) 

Here ni mp , AE± are the actual parameters of the system, and no = 10 12 cm -2 , AE±$ = 4.1 
meV I^Of are the model parameters in our calculations. 



The simplest criterion of the QPC performance is given by average deviation (see (|34])) 

(g n ) = K'jnL ~ 5 ■ 10~ 5 nL( nm). (54) 

For the channel length L = 600(800, 1200) nm and for a criterion of total quantization 
destruction (g) = 0.5, this gives n max w 17(13,8), respectively. The analysis of conduc- 
tance curves shows that the quantization is really destroyed at lower values of (g) ~ 0.2, 
i.e. n max ~ 7(5, 3), which is close to what we see in Fig.|3|. 

The transport and localization lengths can be estimated from ([54]) as follows: 

kr ~ (Kjny 1 ; l loc > l tr ■ n. 
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This leads to the following results: for n open modes k r ~ 20000/n(nm); l\ oc > 20000 (nm). 
These values are much larger than experimentally reasonable contact lengths; they justify 
our theoretical approach of Sec.|TJ. 

In the recent experiment |H] the quantized conductance was observed in an InAs/AlSb 
ballistic constriction with channel length 1000 nm. Up to eight conductance steps were 
detected. Koester et al. suggest that much better performance of their device compared 
to ones utilizing GaAs/ALjGai-^As heterostructures is due to higher ratio of interlevel 
energy spacing to the amplitude of the impurity potential fluctuations. This is in a 
complete agreement with our theory, since the density of quasilocalized states (and the 
parameter k of the theory, see Eq.(f43"D) exponentially drops with growing AE^/a. 

VII. CONCLUSION 

In conclusion, we have investigated the effects of electronic scattering by the soft 
impurity potential in quantum point contacts both with use of an analytical model and 
by direct numerical calculations. 

We have shown that the decisive mechanism of conductance quantization breakdown 
is due to the indirect backscattering of carriers via quasilocalized states at the Fermi 
level. For the realistic contact lengths they can be described in terms of independent 
scatterers, though the electron propagation is coherent. This is due to the smoothness of 
the scattering potentials, leading to very large scattering and localization lengths. 

The performance of the quantum contact is shown to be strictly dependent on the 
ratio of the intermode distance to the amplitude of the random impurity potential (the 
larger the better). 

For the first time we have obtained analytical and empirical formulas both for the 
conductance deviation due to this process and for the probability distribution of these 
deviations in an ensemble of contacts with realistic potentials. The latter has proved to 
be a generalized Poisson distribution. 

The parameters of the distribution obtained numerically agree quite well with ana- 
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lytical calculations based on general assumptions, thus confirming their applicability to 
the case of quantum transport through the QPC in the presence of the random impurity 
potential. 

ACKNOWLEDGMENTS 

We are grateful to Mats Jonson and Ilya Krive for fruitful discussions. 



24 



REFERENCES 



[1] B.J. van Wees, H. van Houten, C.W.J. Beenakker, J.G. Williamson, L.P. Kouwen- 
hoven, D. van der Marel, and C.T. Foxon. Phys.Rev.Lett. 60, 848 (1988); D.A. 
Wharam, M. Pepper, H. Ahmed, J.E.F. Frost, D.G. Hasko, D.C. Peacock, D.A. 
Ritchie, and G.A.C. Jones. J.Phys.C: Solid State Phys., 21, L209 (1988). 

[2] S. Washburn and R.A. Webb, Repts. on Prog, in Physics 55, 1311 (1992). 

[3] L.I. Glazman, G.B. Lesovik, D.E. Khmel'nitskii, and R.I. Shekhter, Pis'ma Zh. Eksp. 
Teor. Fiz. 48, 218 (1988) [JETP Lett. 48, 238 (1991)]. 

[4] J.A. Nixon and J.H. Davies, Phys. Rev. B 41, 7929 (1990); M.J. Laughton, J.R. 
Barker, J.A. Nixon, and J.H. Davies, Phys.Rev. B 44, 1150 (1991); J.A. Nixon, J.H. 
Davies, and H.U. Baranger, Phys.Rev. B 43, 12638 (1991). 

[5] M. Buttiker, Phys.Rev. B 41, 7906 (1990). 

[6] A. Brataas and K.A. Chao, Mod. Phys. Lett. B, 7, 1021 (1993). 

[7] A.M. Zagoskin, R.I. Shekhter, to be published. 

[8] Song He and S. Das Sarma, Phys.Rev. B 48, 4629 (1993). 

[9] D.L. Maslov, C. Barnes, and G. Kirczenow, Phys.Rev. Letters 70, 1984 (1993). 

[10] L.I. Glazman and M. Jonson, Phys.Rev. B 44, 3810 (1991). 

[11] They are only quasilo calized, since their total energy can be equal or exceed the 
energy of propagating states; see below. 

[12] S.A. Gurvitz and Y.B. Levinson, Phys. Rev. B 47, 10 578 (1993). 

[13] J. Imry, in Directions in Condensed Matter Physics, ed. by G.Grinstein and 
G.Mazenko; World Scientific Publ., Singapore, 1986, p. 101. 

[14] T.-Y. Wu and T. Omura, Quantum Theory of Scattering, Prentice-Hall Inc., NJ, 
1962. 

25 



[15] I.M. Lifshitz, S.A. Gredeskul, and L.A. Pastur, Introduction in the theory of disordered 
systems, NY: Wiley, 1988. 

[16] J. Ziman, Models of Disorder, Cambridge University Press, 1979. 

[17] N.K. Patel, J.T. Nicholls, L. Martin- Moreno, M. Pepper, J.E.F. Frost, D.A. Ritchie, 
and G.A.C. Jones, Phys. Rev. B 44, 13549 (1992). 

[18] M.J. Buckingham, Noise in Electronic Devices and Systems, Ellis Horwood Ltd., New 
York- Brisbane- Chichester- Ontario, 1983. 

[19] A.M.Zagoskin and I.O.Kulik, Sov. J. Low Temp. Phys. 16, 533 (1990); E.N.Bogachek, 
A.M.Zagoskin, and I.O.Kulik, Sov. J. Low Temp. Phys. 16, 796 (1990). 

[20] In fact, the value of AE± = 4.1 meV coincides with directly measured in |17|] distance 
between the first and second transverse modes, which ranged from 4.0 to 4.6 meV for 
different gate voltages. 

[21] S. J. Koester, C. R. Bolognesi, E. L. Hu, H. Kroemer, and M. J. Rooks, Phys. Rev. 
B 49, 8514 (1994). 



26 



TABLES 



L(nm)\n 


2 


3 


4 


5 


6 


m 


600 


0.019 


0.018 


0.018 


0.017 


0.016 


0.018 


800 


0.016 


0.017 


0.017 


0.016 


0.015 


0.016 


1200 


0.014 


0.014 


0.014 


0.013 


0.012 


0.014 


7(n) 


0.016 


0.016 


0.016 


0.015 


0.014 


< 7 >= 0.016 



TABLE I. Parameter 7 as determined by the least square fitting of Eq.(51) to the numerical 
data for different contact lengths L and step numbers n. Parameter k is kept equal to (1/320) 
nm _1 for all L, n. 
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FIGURES 

FIG. 1. The model of a quantum point contact. The impurity potential in 2DEG layer. 
Impurity density is 10 12 cm -2 . Isolines are drawn through 1.3 x 1CP 3 eV. The white line shows 
the gate equipotential Ug°\x,y) = Ep (unperturbed contact shape) for L = 800nm. 

FIG. 2. Classification of electronic states in a QPC. E m ±(x) - effective potential in m-th 
mode (ID subband). (a = —2) - reflected states incident from the left; (a = —1) - propagating 
states incident from the left; (a = 0) - quasilocalized states; (a = 1) - propagating states incident 
from the right; (a = 2) - reflected states incident from the right. 

FIG. 3. Conductance vs. gate voltage curves for different realizations and contact lengths 
(see text). Parameter Z = /cf&(0)/7t. a) - contact length 1200 nm, b) - 800 nm, c) - 600 nm. 
Inset: Quantization breakdown parameter g; it equals to the ratio of the area of dashed part of 
the conductance step to its total area. 

FIG. 4. Conductance vs. gate voltage curves for three different realizations for the contact 
length L=10 000 nm (see text). The random potential has been generated on the rectangle 
1290x12900 nm 2 for these calculations. 

FIG. 5. Statistics of conductance fluctuations in QPCs (see text), a) - contact length 1200 
nm, b) - 800 nm, c) - 600 nm. 
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